Energy loss in high energy heavy ion collisions from the Hydro+Jet model 
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We investigate the effect of energy loss of jets in high energy heavy ion collisions by using a 
full three-dimensional space-time evolution of a fluid combined with (mini-)jets that are explicitly 
evolved in space-time. In order to fit the n data for the Au+Au collisions at ^/snn = 130 GeV, the 
space-time averaged energy loss dE/dx{r < 3 fm/c) = 0.36 GeV/fm is extracted within the model. 
It is found that most energy loss occurs at the very early time less than 2 fm/c in the QGP phase 
and that energy loss in the mixed phase is negligible within our parameterization for jet energy loss. 
This is a consequence of strong expansion of the system. 

PACS numbers: 24.85.+p,25.75.-q, 24.10.Nz 



Measurements of high hadrons at the Relativistic 
Heavy Ion Collider (RHIC) and Large Hadron Collider 
(LHC) may provide insight into the excited partonic mat- 
ter, often called a quark gluon plasma (QGP), produced 
in heavy ion collisions [Q. Because jets have to traverse 
the excited matter, their spectra should be changed com- 
pared to the elementary hadron-hadron data. This en- 
ergy loss of jets has been proposed as a possible signa- 
ture of the QGP phase Over the past year, a lot of 
work have been devoted to study the propagation of jets 
through QCD matter |, |, |, § . 

Recently, hadronic transverse momentum distribution 
in Au + Au collisions at y^NN — 130 GeV has been 
measured at RHIC and found that these spectra show 
the depletion at high transverse momentum JtJ. Com- 
parison of the data with pQCD parton model based cal- 
culations 1, | and a transport theoretical study |Tl|l 
shows the indication of energy loss in central Au + Au 
collisions at RHIC energies. Phenomenological stud- 
ies H 12 suggest that the space-time averaged energy 
loss yields small values of dE/dx ~ 0.25-0.3 GeV/fm. In 
addition, jet quenching is found to be the main contribu- 
tion on the elliptic flow parameter in the large transverse 
momentum region Jl3|, Q . Jet energy loss is modeled by 
modifying the fragmentation functions in pQCD parton 
model and usually neglecting the dynamical effects. 

It might be important to consider dynamical effects 
on the jet interactions. The effects of the expansion 
of the system on the energy loss and the azimuthal 
asymmetry have been studied |l^, [l5| [l6|. Hydrody- 
namics provides the space-time evolution of thermalized 
partonic /hadronic matter produced in heavy-ion colli- 
sions g7f |lJl fjl| Recently, one of the authors (T.H.) 
developed a fully three-dimensional hydrodynamic model 
in Cartesian coordinate |2lJ and also in Bjorken coordi- 
nate [jl9| and, moreover, took into account the picture of 
the early chemical freeze-out in this model p2j. Hydro- 
dynamic model calculations at the RHIC energies sug- 
gest that the lifetime of the pure QGP phase is to be 
about 4 fm/c and that the mixed phase exists up to 9-10 



fm/c ]22| |. The density of the system drops rapidly due 
to the longitudinal expansion. 

In this Letter, we explore the dynamical effects on the 
energy loss of jets by taking into account the full 3D 
space-time evolution of a fluid. The hydrodynamic model 
is combined with the jets (Hydro + Jet model) which are 
calculated from the pQCD parton model and explicitly 
propagated in space-time with fluid elements. 

We use full 3D hydrodynamics jl9| in which initial pa- 
rameters are fixed for Au+Au collisions at ^/snn = 130 
GeV [22]| . The transverse profile of the initial energy 
density is assumed to scale with the number of binary 
collisions |fl7f . For 5% central collisions, we choose the 
impact parameter as b = 2.4 fm and the maximum initial 
energy density at the initial time To = 0.6 fm/c as 33.7 
GeV/fm 3 . These parameters lead us to reproduce trans- 
verse momentum spectra of charged hadrons [ p3| up to 
1.0-1.5 GeV/c. It is found that thepT slope is insensitive 
to the thermal freeze-out temperature T when we take 
into account the early chemical freeze-out 22 . We fix 
T th = 140 MeV throughout this Letter. 

We include hard partons using pQCD parton model, 
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where Y\ and Y 2 are the rapidities of the scattered par- 
tons and x\ and x 2 are the fractions of momentum of 
the initial partons. The parton distribution functions 
f a (x, Q 2 ) are taken to be CTEQ5 leading order ||. We 
use Q 2 = p T for the evaluation of parton distribution. 
The minimum momentum transfer pT.min = 2.0 GeV/c 
is assumed. The summation runs over all parton species 
and relevant leading order QCD processes 



(2) 
(3) 
(4) 

are included in addition with the initial and final state 
radiation to simulate the emission of multiple soft gluons. 



q + q -> q + q , q + q^ q + q 
q + q ->■ g + g q + g^ q + g, 
g + g -> q + q, g + g -> g + g. 
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FIG. 1: Comparison of tt spectra from hydro + jet with the 
energy loss parameter ea = 0.06 GeV fm 2 (solid line), hydro 
+ jet with ea = 0.3 GeV fm 2 (dashed line), and hydro only 
(dotted line) for b = 3.35 fm to the PHENIX data for 
centrality of 10%. The result of hydro + jet for dE/dx = 0.2 
GeV/fm calculation is shown in open squares. 

Gaussian primordial transverse momentum /ct distribu- 
tion with the width of (k T ) = 1 GeV 2 is assigned to 
the shower initiator in the QCD hard 2-^2 processes. 
We use PYTHIA 6.2 [p5j to simulate each hard scatter- 
ing in the actual calculation. A factor K is used for the 
higher order corrections. We chose K = 2 to fit the UA1 
data of pp at yfs = 200 GeV [^6) . In order to convert 
hard partons into hadrons, we use an independent frag- 
mentation model using PYTHIA after hydro simulations. 
We have checked that this hadronization model provides 
good agreement with the transverse spectra of charged 
hadrons above pr = 1 GeV/c in pp data |26| . 

The number of hard partons is assumed to scale with 
the number of hard scattering which is estimated by 
using Woods-Saxon nuclear density. This assumption 
is consistent with the peripheral Au+Au collision at 
\J snn = 130 GeV |J. The space coordinate of a parton 
in longitudinal direction is taken to be rj s = Y, where 
r/ s = (l/2)log[(t + z)/(t — z)], and Y is a momentum 
rapidity. The transverse coordinate of a parton is spec- 
ified by the number of binary collision distribution for 
two Woods-Saxon distributions. Most hard partons are 
produced with the formation time of ~ 1/.PT assuming 
the uncertainty relation. Since the initial time of the hy- 
drodynamical simulation is tq = 0.6 fm/c, partons are 
assumed to travel freely up to r . We neglect the nu- 
clear shadowing effect in the present work for simplicity, 
because it is small at high transverse momentum [g, [12| . 

We assume that the form of energy loss is simply 

dE 

— = eap(T,r), (5) 

where e is an energy loss per scattering and a is a 
parton-parton cross section. From a hydrodynamic sim- 



FIG. 2: Comparison of ir° spectra from hydro + jet (solid line) 
and hydro only (dotted line) with the energy loss parameter 
ea = 0.06 GeV fm 2 and b = 12.1 fm to the PHENIX data §] 
for peripheral collision (60-80%). 

ulation, we obtain the space-time evolution of temper- 
ature T(t, r). A parton density p(r,r) in the QGP 
phase is calculated from T(r, r). For the mixed phase 
at T = T c (=170 MeV), we put p(r,r) = / qgp (t, r)p(T c ) 
where the fraction of the QGP phase in the mixed phase 
fQGp(T 7 r) is calculated from energy density e(r, r) and 
the maximum and minimum energy densities in the 
mixed phase £qgp and £had as 

, , > e(r, r) - e had 
/qgp (T,r) = . (6) 

£QGP — ^had 

We assume that the energy loss for a quark jet is half 
of the energy loss of a gluon p7j . In the present work, 
we neglect the possible space-time variation of parton- 
parton cross section and regard the product ea as an 
adjustable free parameter in the model. Feedback of the 
energy to fluid elements is ignored, because we use the 
relatively small values of energy loss parameter used in 
the work. 

In Fig. |l|, the transverse momentum distributions of 
neutral pions for central 10% Au + Au collisions at 
y/IWN = 130 GeV from PHENIX experiment § are com- 
pared to the results of the hydro + jet model with dif- 
ferent energy loss parameters ea — 0, 0.06 and 0.3 GeV 
fm 2 together with the hydro result by dotted line. In the 
calculation, we choose the impact parameter b — 3.35 
fm for 0-10% central events. The corresponding num- 
ber of binary collisions is Aginary = 906. In order to 
connect the contribution from jets and a fluid smoothly, 
low pt pions from jets are cut by a switch function 
{1 + tanh[3(p T - 1.5GeV)]}/2 @. We reproduce p T 
spectrum for ir° by choosing ea — 0.06 GeV fm 2 , while 
the result which is not taken into account a energy loss 
overestimates the neutral pion spectrum. The hydro + 
jet result for dE/dx — 0.2 GeV/fm which is obtained 
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FIG. 3: The ratio Raa for neutral pions in Au+Au collisions 
with b = 3.35 fm. The error bars of the PHENIX data jjj] 
contain the statistical error only. 



by neglecting the density dependence in Eq. (||) is also 
shown in Fig. [|. We note that the hydro + jet result for 
en- = 0.3 GeV fm 2 is reproduced by the constant energy 
loss of dE/dx = 1.0 GeV/fm. As shown in Fig. |, tt 
spectra for peripheral (60-80%) Au + Au collisions 
can be reproduced only by changing the impact param- 
eter to b = 12.1 fm (Aginary = 20) and leaving the other 
parameters. The hydro -I- jet result in peripheral colli- 
sions does not change when energy loss is not included, 
because the parton density and the volume of the QGP 
phase are small. On the other hand, the hydro + jet 
model overpredicts the data without energy loss in the 
central collisions. The deviation from data in the trans- 
verse momentum range below 1.5 GeV/c in the peripheral 
collisions might be reasonable, because we do not expect 
that local thermalization is achieved in such a peripheral 
collision and that hydrodynamics is reliable. 

The jet quenching is quantified by the ratio of the par- 
ticle yield in A+A collisions to the one in p+p collisions 
scaled up by the number of binary collisions 



Raa = 



d 2 N A+A /dp T dr] 
NbuuuydPNP+P/dprdT)' 



(7) 



The transverse momentum dependence of Raa can pro- 
vide information about the mechanism of jet quenching 
p8[ . In Fig. H the result from the Hydro+Jet model with 
ecr = 0.06 GeV fm 2 is compared to the PHENIX data @. 
Our result is almost flat in the range 2 < px < 4 GeV/c. 

Let us now turn to the study of the dynamical effects 
on the jet energy loss. The maximum and the average 
thermalized parton densities at 77 s = from hydrody- 
namic simulations are plotted as a function of proper 
time in the collision of impact parameter b — 3.35 (12.1) 
fm in upper (lower) panel of Fig. |^. As seen in the fig- 
ures, thermalized parton density drops rapidly due to the 
strong longitudinal expansion of the system produced in 
heavy ion collisions. We expect that jets are likely to loss 
their energies in the time span of less than r ~ 2 fm/c. 
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FIG. 4: Maximum (circles) and average parton density 
(squares) as a function of proper time at impact parameter 
b = 3.35 fm (upper) and b — 12.1 fm (lower). 



In order to see when jets lose their energies, we plot in 
Fig. ||, the ratio of the numbers of jets N(t)/N(t ) at 
Pt = 5 GeV/c. It is found that about 90% of jet quench- 
ing occurs up to r ~ 2 fm/c, although the QGP phase 
exists up to t ~ 4 fm/c and the mixed phase lasts until 
t ~ 10 fm/c within the hydro parameters used in this 
work. This holds true, even when we take the total den- 
sity neglecting the existence of the hadrons in the mixed 
phase in Eq. (^) in order to get the maximum energy 
loss. Therefore, within our model, only QGP phase is 
responsible for the jet energy loss. 

We showed in Fig. [I] the energy loss dE/dx — 0.2 
GeV/fm gives the same amount of energy loss as the 
energy loss parameter ecr = 0.06 GeV fm 2 . Let us es- 
timate a space-time averaged energy loss parameter for 
central collision (b = 3.35 fm). If we take a space-time 
average parton density p(r), p(r < 1.0 fm/c) = 11 fm , 
p{r < 2.0 fm/c) = 7.7 far 3 , and p(r < 3.0 fm/c) = 6.0 
fm -3 , then the space-time averaged energy loss yields 
dE/dx(r < 1.0 fm/c) ~ 0.66 GeV/fm, dE/dx(r < 
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FIG. 5: Jet quenching rate defined as the ratio of the numbers 
of jets N(t)/N(to) at pt = 5 GeV/c as a function of proper 
time for b = 3.35 fm. The energy loss parameter ea = 0.06 
GeV fm 2 is used. 



2.0 fm/c) ~ 0.46 GeV/fm, and dE/dx(r < 3.0 fm/c) ~ 
0.36 GeV/fm with the best fit parameter of ea = 0.06 
GeV fm 2 . Those values are close to that of the result 
calculated by the pQCD parton model ||, in which the 
effect of expansion is not included. The reason is that 
the energy loss parameter can be regarded as an average 
of roughly a short time of r < 2 fm/c. 

In summary we proposed the full 3D hydrodynamical 
model combined with (mini-) jets, where jets are explic- 
itly propagated in space-time with the hydro simulation. 
In particular, this model allows us to study the dynamical 
effects on the jet energy loss. We estimated the energy 
loss of partons by the fluid elements whose temperature 



is above the critical value T c . It is found that energy loss 
of jets occurs in the pure QGP phase of r < 2 fm/c and 
that the contribution of energy loss in the mixed phase 
is negligible under the assumption of the energy loss for- 
mula Eq. (|^). This indicates that suppression of hadronic 
high px spectra contain information about the early stage 
of partonic matter. However, it has been shown that the 
formula for the energy loss has non-trivial energy depen- 
dence H [| and it is found that energy loss is sensitive to 
the critical point [^9). It is very important to take into 
account the coherent (Landau-Pomeranchuk-Migdal) ef- 
fect in the calculation of radiation spectrum. A result 
obtained by using a different energy loss formula will be 
presented elsewhere. 

A prediction at higher transverse momentum region 
is under progress. Elliptic flow parameter should be 
also studied within the model. It should be studied 
the charged particle spectra in order to have the uni- 
fied understanding of the jet quenching mechanism in 
the medium. 

In this study, we do not include the energy loss effects 
before hydrodynamical evolution. It might be interesting 
to study to what extent partons lose energy before ther- 
malization using non-equilibrium models [fill |30(| . Be- 
cause parton density is maximum in this stage, energy 
loss by these partons should be important. 
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